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Introduction 


The  long-term  goal  of  this  research  projec  t  is  to  develop  CCR5  antagonists  with  the 
structure  feature  of  anibam  ine  as  novel  anti  pr  ostate  cancer  agents  .  T  he  tasks  of  this 
project  include:  1.  Chemical  synthesis  of  the  ligand  designed  as  anibamine  derivatives;  2. 
Radio-ligand  competition  binding  assay  to  ev  aluate  the  b  inding  affinity  of  the  ligands 
synthesized;  3.  Characterization  of  CCL5  and  CCR5  expre  ssion  in  SV40T-immortalized 
human  prostate  epithe  lial  ce  11  lines  of  the  sam  e  cellular  lin  eage  with  v  arying 
tumorigenicity  and  m  etastatic  capacity  in  vivo  (P69,  M2182,  M2205  and  M12);  4. 
Investigation  of  the  impact  of  chemokine  receptor  CCR5  antagonist  on  prostate  cancer 
cell  growth  and  progression  using  M  12,  PC-3,  DU-145  and  LNCa  P  cell  lines;  5. 
Molecular  modeling  study.  In  order  to  achiev  e  these  aim  s,  we  have  first  designed  and 
synthesized  a  series  of  anibam  ine  anal  ogs  based  on  “Deconstruction-Reconstruction- 
Elaboration”  method.  So  far  thirty-two  novel  ligands  have  been  prep  ared  as  derivatives 
of  aniba  mine  for  biological  screening  assa  ys.  Second,  several  progressive  biological 
assays  to  ev  aluate  the  anti  pros  tate  cancer  activity  of  the  ligands  have  been  set  up  and 
pursued.  Prim  arily,  a  RANTES-b  inding  inhibition  assay  protoc  ol  has  been  set  up  and 
ready  to  te  st  our  drug  candidates.  Additionally,  biological  s  creening  includes  the  te  st  of 
capacity  of  these  novel  compounds  to  inhibit  proliferation  and/or  apoptosis  by  the  human 
prostate  cancer  cell  lines  M12,  PC-3,  DU-145  and  LNCaP  has  been  conducted  to  evaluate 
their  efficacy.  Also  we  have  assessed  cha  nges  in  the  expression  of  CCL5  and  CCR5  in 
human  prostate  epithelial  cell  lines  of  the  same  cellular  lineage  but  with  differing  in  vivo 
phenotypes  (P69SV40TAg,  M2 182,  M2205,  and  M  12).  A  molecular  m  ode  ling  study  has 
been  conducted  to  elucidate  the  interaction  of  anibam  ine,  the  lead  com  pound,  to  the 
target  protein,  chem  okine  receptor  CCR5,  to  further  facilitate  our  next  generation 
molecular  design.  This  part  of  work  has  been  published  in  the  Journal  of  Che  mical 

Informatics  and  Modeling. 
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Annual  Progress  Report 

The  long-term  goal  of  this  research  project  is  to  develop  CCR5  antagonist  with  the  structure  feature  of 
anibamine  as  novel  anti  prostate  cancer  agents.  Based  on  the  approved  Statement  of  Work  of  this  award,  the 
following  tasks  will  be  accomplished  within  two-year  period  of  time: 

Task  1.  Chemical  synthesis  of  the  ligand  designed  as  anibamine  derivatives 

Task  2.  Radio-ligand  competition  binding  assay  to  evaluate  the  binding  affinity  of  the  ligands  synthesized 

Task  3.  Characterization  of  CCL5  and  CCR5  expression  in  SV40T-immortalized  human  prostate  epithelial  cell 
lines  of  the  same  cellular  lineage  with  varying  tumorigenicity  and  metastatic  capacity  in  vivo 

Task  4.  Investigation  of  the  impact  of  chemokine  receptor  CCR5  antagonist  on  prostate  cancer  cell  growth  and 
progression  using  M12,  PC-3,  DU-145  and  LNCaP  cell  lines 

Task  5.  Molecular  modeling  study 

In  the  past  one  year  all  five  tasks  have  pursued  accordingly  and  the  results  are  summarized  below: 


Task  1.  Chemical  synthesis  of  the  ligand  designed  as  anibamine  derivatives 


Our  laboratory  has  accomplished  the  total  synthesis  of  anibamine  and  published  the  results  in  Organic  Letters, 
2007,  9(10),  2043-6.  This  biosynthetic  approach  is  being  applied  directly  to  our  current  chemical  synthesis  of 
the  newly  designed  derivatives  of  anibamine.  The  initial  structural  modification  of  anibamine  is  following  the 
“Deconstruction-Reconstruction-Elaboration”  method,  as  proposed  in  Figure  1. 


Decontruction:remove 
this  side  chain  to  test 
its  influence  on  the 
anti-cancer  activity 


Deconstruction:  remove 
the  ring  system  to  see  its 
influence  on  activity 


Elaboration:  different  configuration:  cis 
or  trans;  or  saturated  side  chain 

Decontruction:remove 
this  side  chain  to  test 
its  influence  on  the 
anti-cancer  activity 

Elaboration:  different  configuration:  cis  or 
trans;  or  saturated  side  chain 


Figure  1.  Initial  structure  modification  of  anibamine 


1.1  Deconstruction  procedure 


1.1.1  Core  ring  system  The  central  ring  system  of  anibamine  largely  determines  the  molecular 
conformation  of  anibamine.  Thus,  by  removing  the  fused  ring,  the  conformation  of  the  molecule  will  be 
changed.  Shown  in  scheme  1,  the  synthesis  of  non-ring  derivatives  of  anibamine  was  very  straightforward.  The 
methyl  group  was  introduced  with  little  difficulty. 


OMe 

Scheme  2.  The  synthesis  of  anibamine  derivative  with  non-fused-ring  system 


Another  obvious  modification  is  to 
remove  the  core  ring  system  totally 
and  test  if  that  will  influence  the 
activity.  The  synthesis  of  the  ligand 
was  accomplished  by  following  the 
routes  in  scheme  2. 
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1.1.2  Side  chains  Alternative  removal  of  side  chain  at  position  5  or  3  will  test  the  influence  of  side  chains 
to  the  anti  prostate  cancer  activity  of  anibamine.  The  synthesis  of  these  two  compounds  was  finished  by 
following  scheme  3  and  4. 


Scheme  3.  Synthesis  of  one  side  chain  at  position-3  derivative  of  anibamine 


H 


nn  ncch2cooh 


h3c^chcn 
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CH3C(OEt)3 
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O  O 


OEt 


NC 


P2O5/TBAB  NC 
Toluene 


N'  OH  N  Br 

Scheme  4.  The  synthesis  of  one  side  chain  at  position  5  derivative  of  anibamine 


1.2  Reconstruction  procedure  Based  on  the  biological  screening  results  (being  conducted  right  now)  of 
the  deconstruction  products  (1.1),  we  will  reconstruct  the  skeleton  of  the  active  molecules  to  verify  the  essential 
pharmacophore(s)  in  the  parent  compounds.  Total  of  six  compounds  may  be  synthesized  to  verify  their  activity. 


1.3  Elaboration  procedure 


1.3.1  Double  bond  configuration  We  further  changed  the  double  bond  configuration  from  cis  to  trans  in  the 
molecule  to  test  such  elaboration’s  influence  on  the  anti  prostate  cancer  activity  (Scheme  5). 


Scheme  5.  Anibamine  derivative  synthesis  with  different  double  bond  configuration 


1.3.2  Core  ring  system 
modification  By  changing  the 
size  of  the  fused  ring,  the 
conformation  of  the  molecule 
will  be  changed.  This  may 
change  its  binding  affinity  to 
CCR5  significantly  and  further 
influence  their  anti  prostate 
cancer  activity.  As  shown  in 
scheme  6,  the  syntheses  of  ring- 
size  modified  derivatives  of 
anibamine  were  straightforward. 


Scheme  6.  Synthesis  of  ring  system  modified  derivatives  of  anibamine 


1.3.3  Side  chain  further  modification  Based  on  the  fact  that,  first,  the  calculated  log  Kow  for  anibamine  is  9.1, 
which  indicates  that  its  lipophilicity  is  significant  higher  than  the  value  set  forth  by  “Lipinski’s  rule  of  5”  for 
drug-like  compounds;  second,  in  comparing  the  chemical  structure  of  anibamine  with  other  known  CCR5 
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antagonists  the  major  difference  is  that  the  two  anibamine  side  chains  are  simple,  undecorated,  aliphatic  chains, 
we  decided  to  introduce  more  chemical  structure  features  on  these  side  chains  that  will  effectively  drive  this 
structure-activity  relationship  study.  We  replaced  the  long  aliphatic  side  chains  with  substituted  benzene  ring 
systems  with  the  same  or  similar  number  of  carbons.  Some  of  the  model  compounds  are  listed  in  Table  1.  The 
chemical  synthesis  of  these  compounds  was  conducted  by  following  scheme  1,  3  and  4. 

Table  1.  Spacer  and  substitutions  for  the  side  chain  modification 

r  Spacer  X  None,  CH2,  CH2CH2 

R  (bulky  effect)  OCH3,  OCH2CH3,  O-i-Pr,  O-c-pentyl,  O-c-hexyl 

R  (electronic  effect)  NH2,  OCH3,  CH3,  H,  F,  Cl,  Br,  N02,  CF3,  S03Me 

Among  all  the  ligands  described  in  Table  1,  we  have  finished  the  synthesis  of  twenty-four  compounds.  We 
hope  to  be  able  to  apply  the  biological  screening  results  of  these  ligands  to  guide  our  choices  in  our  side  chain 
modification  (up  to  forty-five  novel  ligands  may  be  synthesized  eventually). 

In  summary,  total  of  thirty-two  novel  ligands  as  anibamine  derivatives  have  designed  and  synthesized  through 
multi-step  chemical  synthesis.  They  all  have  been  fully  characterized  by  NMR,  IR,  and  MS.  Elementary 
analysis  were  also  conducted  for  each  compound  for  final  verification  of  their  chemical  composition. 

Task  2.  Radio-ligand  competition  binding  assay  to  evaluate  the  binding  affinity  of  the  ligands  synthesized 

2.1  RANTES  binding  inhibition  assay  We  have  setup  a  radio-ligand  binding  assay  protocol  to  conduct  the 
RANTES  binding  inhibition  assays  to  test  the  CCR5  receptor  affinity  of  all  the  compounds  synthesized.  The 
assay  is  now  being  conducted  for  all  the  compounds  so  far  synthesized. 

To  test  the  CCR5  receptor  affinity  of  the  compounds  synthesized,  Chinese  hamster  ovary  (CHO)  cells 
expressing  human  CCR5  will  be  used  in  the  binding  experiments.  The  experiment  will  be  conducted  at  room 
temperature  in  the  presence  of  [125I1RANTES.  IC50  values  (the  concentration  required  to  inhibit  the  binding  of 
[  IJRANTES  by  50%)  will  be  used  to  evaluate  the  relative  affinity  of  the  compounds  to  CCR5.  Normally, 
compounds  with  an  IC5o  concentration  at  lOpM  or  lower  are  considered  active  and  will  be  further  evaluated  in 
additional  in  vitro  assays.  The  binding  assay  will  provide  preliminary  data  for  further  molecular  design, 
especially  as  guided  by  molecular  modeling  where  downstream  biological  effects  are  not  explicitly  considered. 

2.2  Protocol 

2.2.1  CCR5/CHO  Cell  Transfection  CHO  Cells  are  transfected  with  CCR5  plasmid  DNA  (UMR  cDNA 
Resource  Center)  using  LipofectAMINE  Reagent  (GibcoBRL)  according  to  the  manufacturer’s 
recommendation.  Briefly,  a  lipid/DNA  mixture  is  prepared  in  200  pL  of  serum-free  media  (3  pL  lipid  and  1  pg 
DNA)  and  incubated  for  45  min  at  22°C.  Serum-free  media  is  then  added  to  bring  the  final  volume  to  1  mL, 
which  is  then  added  to  each  well  of  cells.  After  3  hr,  the  lipid/DNA  mixture  is  removed  and  replaced  with  1  mL 
of  FBS-containing  medium.  For  selection  of  stable  cell  lines,  the  media  is  replaced  with  selection  media 
(containing  4  to  6  mg/mL  of  G418)  after  36  hr.  These  cells  are  then  split  into  96  well  plates  at  a  concentration 
of  1  cell  per  three  wells.  The  entire  time  the  cell  are  being  grown  in  the  presence  of  1  mg/mL  of  G418  to 
continue  to  maintain  the  selection  process.  Once  the  cells  get  confluent,  Western  blots  will  be  conducted  for 
each  clone  utilizing  a  human  CCR5  polyclonal  antibody.  The  cloned  cells  with  CCR5  expression  will  be 
adopted  to  further  cell  culture. 

2.2.2  Cell  Culture  and  Drug  Treatments  Cells  are  cultured  in  a  humidified  atmosphere  of  5%  C02/95%  air 
at  37°C.  The  medium  is  a  1:1  mixture  of  DMEM  and  Nutrient  Mixture  F12  containing  5%  fetal  bovine  serum, 
100  units/mL  each  of  penicillin  and  streptomycin.  Transfected  cell  medium  also  contains  0.4  mg/mL  G-418. 

2.2.3  Membrane  Homogenate  Preparation  Cells  are  harvested  by  replacing  the  media  with  PBS  +  0.4% 
EDTA  and  collected  by  centrifugation  at  1000  x  g.  Samples  are  homogenized  in  20  volumes  ice-cold 
membrane  buffer  (50  mM  Tris-HCl,  3  mM  MgCl2j  1  mM  EGTA,  pH  7.4)  with  a  Polytron.  The  homogenate  is 
centrifuged  at  48,000  x  g  for  10  min,  resuspended  in  assay  buffer  (50  mM  Tris-HCl,  3  mM  MgCF,  0.2  mM 
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EGTA,  pH  7.4),  centrifuged  again  as  above,  and  resuspended  in  assay  buffer.  Protein  values  are  determined  by 
the  Bradford  method. 

2.2.4  CCR5  Receptor  Binding  Saturation  binding  is  performed  by  incubating  membranes  for  90  minutes 

at  30°C  with  0.5-15  nM  [125I]RANTES  in  assay  buffer  in  a  0.5  mL  volume.  Non-specific  binding  is  determined 
with  200  nM  RANTES.  For  competition  assays,  membranes  are  incubated  as  above  with  2  nM  [  IJRANTES 
and  various  concentrations  of  unlabeled  ligand,  to  determine  competitor  affinity  for  CCR5.  The  reaction  is 
terminated  by  rapid  filtration  (GF-C  filter  paper  is  presoaked  for  at  least  half  hour  in  0.3%  PEI  solution). 
Bound  radioactivity  is  determined  by  gamma  counter  for  [125I]  immediately  [54], 

2.2.5  Data  analysis  For  competition  binding  assay,  linear  regression  analysis  of  Hill  plots  (plot  of  log[B  /  (Bo 
-  B)]  as  a  function  of  log[X];  where  B  is  the  pmol/mg  of  radioligand  bound  and  X  is  the  concentration  of 
unlabeled  competitor)  will  be  conducted  to  determine  the  IC50  value  for  each  ligand.  “B”  will  be  calculated  as 
%  of  radioligand  bound  at  each  concentration  of  competitor  and  Bo  will  be  normalized  to  be  100%.  The  IC50 
values  (inverse  log  X  where  B  =  50%)  will  then  be  determined  and  corrected  to  K,  values  using  the  Cheng- 
Prusoff  equation:  K;  =  IC50  /  (1  +  (L/KD)),  where  L  is  the  concentration  and  KD  is  the  KD  value  of  the 
radioligand.  All  linear  and  non-linear  curve-fitting  analyses  will  be  performed  using  Prism  4.0  software. 

2.3  Results  The  binding  assay  is  now  being  conducted  under  the  supervision  of  Dr.  Dana  E.  Selley  of  the 
Department  of  Pharmacology  and  Toxicology  at  Virginia  Commonwealth  University.  Both  the  Pi’s  lab  and  Dr. 
Selley ’s  lab  have  fully  functional  binding  assay  facility.  The  assays  will  be  pursued  mainly  in  the  Pi’s  lab  while 
Dr.  Selley’ s  lab  will  provide  necessary  technical  support.  The  results  will  be  summarized  and  reported  in  due 
course. 

Task  3.  C  haracterization  of  CCL5  and  CCR5  exp  ression  in  S  V40T-immortalized  human  prostate 
epithelial  cell  lines  of  the  same  c  ellular  lineage  with  varying  tumo  rigenicity  and  metastatic  capacity  in 
vivo 

SYBR-Based  Real-time  PCR:  protocol  and  results 

Gene  expression  of  CCR5  and  CCL5  was  quantitated  by  SYBR-based  Real-time  PCR.  The  U6  gene  was  used 
as  internal  control.  cDNA  was  synthesized  using  the  iScript  cDNA  Synthesis  Kit  (Biorad)  in  a  20  pi  reaction 
volume  including  4  pi  of  iScript  reaction  mix,  1  pi  of  iScript  reverse  transcriptase  and  1  pg  of  total  RNA.  The 
reactions  were  incubated  first  at  25°C  for  5  min  and  then  at  42°C  for  30  min  followed  by  inactivation  at  85°C 
for  5  min.  The  CCL5,  CCR5  and  U6  genes  were  amplified  by  Real-time  PCR  in  a  7900HT  Sequence  Detection 
System  (Applied  Biosystems).  The  20-pl  PCR  mixture  included  10  pi  of  the  2X  Universal  SYBR  Master  Mix 
(Applied  Biosystems),  1  pi  of  primers  (10  pM)  and  3  pi  of  reverse  transcription  (RT)  product.  The  reactions 
were  incubated  in  a  96-well  plate  at  50°C  for  2  min,  95°C  for  10  min,  followed  by  40  cycles  of  95°C  for  15  s, 
55°C  for  35  s  and  68°C  for  35  s.  Primers  used  for  the  CCL5  gene  were  5’- 
CTC  ATTGCTACTGCCCTCTGCGCTCCTGC-3  ’  and  5  ’  -GCTC  ATCTCC  AAAG  AGTTGATGTACTC-3  ’ . 
Primers  used  for  CCR5  gene  were  5’-TTGAGTCCGTGTCACAAGCCC-3’  AND  5’- 
AATAATTGCAGTAGCTCTAACAGG-3’  Primers  used  for  U6  gene  were  5’-CTCGCTTCGGCAGCACA-3’ 
and  5’-AACGCTTCACGAATTTGCGT-3’.  The  threshold  cycle  (CT)  was  defined  as  the  fractional  cycle 
number  at  which  the  fluorescence  passed  the  fixed  threshold  (0.2).  The  relative  quantity  (RQ)  of  the  target 
miRNAs  was  estimated  by  the  AACT  study  using  U6  expression  as  an  endogenous  control  for  each  reaction. 
Relative  expression  was  calculated  using  the  comparative  CT  method. 
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CCL5 


CCR5 


Figure  2.  The  qRT-PCR  results  to  characterize  the  expression  of  CCL5  and  CCR5  in  different  cell  lines. 

From  Figure  2,  it  appears  that  the  expression  levels  of  CCL5  and  CCR5  is  correlated  with  the  tumorigenicity 
and  metastatic  capacity  of  different  cell  lines.  That  is,  the  higher  the  tumorigenicity  and  metatstatic  capacity  the 
cell  line  carries,  the  higher  the  expression  level  of  its  chemokine  receptor  CCR5  would  be.  Further 
interpretation  and  verification  of  the  results  will  be  conducted  after  more  prostate  cancer  cell  lines  would  be 
examined. 

Task  4.  Investigation  of  the  impact  of  chemokine  receptor  CCR5  ant  agonist  on  prostate  cancer  cell 
growth  and  progression  using  M12,  PC-3,  DU-145  and  LNCaP  cell  lines 

4.1  Anti  Proliferation  Assays  The  capacity  of  each  synthesized  compounds  to  inhibit  the  proliferation  of 

several  prostate  cancer  cell  lines  will  be  tested  primarily  as  the  first  biological  assay  to  evaluate  their  anti  cancer 
activity. 

4.2  Protocol 

All  cell  lines,  P69,  PC-3,  DU-145  and  M12,  were  incubated  at  37  °C  in  the  presence  of  5%  C02.  RPMI  1640 
serum  free  media  (GIBCO  Invitrogen)  containing  1  %  L-glutamine,  0.1  %  ITS  (insulin,  5pg/ml;  transferrin, 
5pg/ml;  and  selenium,  5  pg/ml;  Collaborative  Research,  Bedford)  and  0.1  %  gentamicin  was  used  to  cultivate 
all  cells. 

The  media  for  DU- 145  and  PC-3  cell  lines  also  included  10  %  fetal  bovine  serum  (FBS).  The  serum  containing 
media  for  the  M12  and  P69  cell  lines  contained  5  %  FBS.  While  the  DU-145,  PC-3  and  P69  cells  were  all 
originally  plated  in  media  containing  serum,  the  M12  cells  were  originally  plated  in  serum  containing  media 
and  then  the  next  day  were  plated  in  serum  free  media  and  0.01  percent  epidermal  growth  factor  (EGF).  B. 
Anti-Proliferation  Drug  Assays 

The  effect  of  compounds  designed  as  CCR5  antagonists  on  various  prostate  cancer  tumor  cell  lines  was 
assessed  utilizing  WST-1  Cell  Proliferation  Reagent  (Roche).  Prostate  cancer  tumor  cells  (DU-145,  PC-3,  M12 
and  P69)  were  plated  out  into  96  well  plates  (BD  Falcon,  VWR)  at  a  concentration  of  2000  cells  per  well.  For 
the  purpose  of  this  study  all  cells  were  plated  in  RPMI  media  (GIBCO  Invitrogen)  containing  10  %  fetal  bovine 
serum,  1  %  L-glutamine,  0.1  %  ITS  and  0.1  %  gentamicin.  The  cells  were  placed  in  a  total  volume  of  100 
microliters  per  well  and  incubated  overnight  at  37  °C  and  5%  C02.  Various  concentrations  of  drug  were  plated 
the  following  day,  maintaining  a  final  total  well  volume  of  150  microliters  per  well.  In  wells  containing  only 
cells  and  media  or  just  media,  a  volume  of  50  milliliters  of  PBS  was  added  to  each  well.  All  drugs  were 
dissolved  in  a  minimum  amount  of  dimethysulfoxide  (Sigma-Aldrich)  and  the  further  dissolution  was 
conducted  with  PBS  buffer  (GIBCO  Invitrogen).  Once  the  drugs  were  plated  they  were  allowed  to  incubate 
with  the  cells  for  24  hours  at  37  °C  and  5%  C02.The  next  day  10  microliters  of  WST-1  Cell  Proliferation 
Reagent  (Roche)  was  added  to  each  well,  and  incubated  at  37  °C  and  5%  C02  for  3  hours.  The  absorbance  of 
each  well  was  then  measured  using  an  EL  312e  Microplate  Bio-kinetics  Reader  (BIO-TEK  Instruments). 
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Calculations  of  %  inhibition  and  IC50  of  each  drag  was  performed  utilizing  a  spreadsheet  (Microsoft  Office 
Excel  2007).  The  percent  inhibition  was  calculated  by  subtracting  the  average  absorbance  of  the  cells  in  the 
presence  of  drag  from  the  absorbance  of  the  cells  with  just  media  and  then  this  value  was  divided  by  the 
difference  between  the  absorbance  of  the  cells  without  the  drag  and  the  absorbance  of  the  media,  all  of  which 
were  multiplied  by  one  hundred.  The  equation  is,  percent  inhibition  =  [  (Acells  -  Adrag)  /  (Acells  -  Amedia)  ]  x 
100,  where  A  is  the  absorbance,  cells  is  wells  growing  cells  with  no  drag,  drag  is  cells  growing  in  the  presence 
of  drag  and  media  is  just  media  in  the  absence  of  both  cells  and  drag. 

4.3  Results  and  Discussion  So  far,  nine  compounds  have  been  evaluated  fully  in  three  prostate  cancer  cell 
lines  (DU- 145,  PC-3,  and  M12)  for  three  timelines  (24  hours,  48  hours,  and  72  hours)  at  a  serious  of 
concentration  (from  1  nM  to  10  uM).  As  indicated  in  Figure  3,  the  inhibition  of  proliferation  effect  varies 
significantly  depending  on  the  chemical  structure  of  the  ligands.  This  provides  important  Stracture-Activity- 
Relationship  information  to  direct  us  for  the  next  step  molecule  design.  Further  evaluation  for  all  other  ligands 
(synthesized  and  being  synthesized)  will  be  conducted  following  the  same  protocol.  Once  a  complete  data 
available,  Quantitative  Structure  Activity  Relationship  (QSAR)  will  be  conducted  to  provide  detailed 
information  for  the  next  step  molecular  design. 
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Figure  3.  Inhibition  effect  of  prostate  cancer  cell  line  proliferation  by  new  ligands  synthesized. 
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Task  5.  Molecular  modeling  study 

5.1.  Homology  modeling  study  of  CCR5  receptor  3D  model 

As  the  first  CCR5  natural  product  ligand,  anibamine  provides  a  novel  structural  skeleton  that  is  remarkably 
different  from  all  lead  compounds  previously  identified  as  CCR5  antagonists.  In  order  to  characterize  its 
interaction  with  the  target  protein,  we  conducted  comparative  docking  studies  of  anibamine  with  several  other 
known  CCR5  antagonists  in  two  CCR5  homology  models  built  based  on  the  crystal  structures  of  bovine 
rhodopsin  and  human  p2-adrenergic  receptor.  From  our  study,  the  binding  pocket  of  anibamine  has  some 
common  features  shared  with  other  high  affinity  CCR5  antagonists,  suggesting  that  they  may  bind  in  similar 
binding  sites  and/or  modes.  At  the  same  time,  several  unique  binding  features  of  anibamine  were  identified  and 
it  will  likely  prove  beneficial  in  future  molecular  design  of  novel  CCR5  antagonists  based  on  the  anibamine 
scaffold.  The  part  of  research  results  have  been  published  in  the  Journal  of  Chemical  Informatics  and 
Modeling,  2009,  49(1),  120-32. 

5.2  Small  molecule  construction  For  all  the  small  molecules  (designed  and  synthesized  ligands),  their 
conformation  has  been  optimized  after  being  constructed  in  Sybyl  7.0.  They  are  ready  for  the  next  step 
molecular  modeling  study. 
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Key  Research  Accomplishments 


•  A  series  of  novel  ligands  have  been  designed  based  on  the  chemical  structure  of 
anibamine,  the  first  nature  product  chemokine  receptor  CCR5  antagonist. 

•  Chemical  synthesis  of  thirty-two  (32)  novel  ligands  among  this  series  has  been 
finished  and  all  the  ligands  have  been  characterized  fully  to  verify  their  chemical 
structures. 

•  A  RANTES  binding  inhibition  assay  has  been  set  up  and  the  protocol  has  been 
verified  using  RANTES  and  TAK-779.  The  method  is  ready  for  the  screening  of 
synthesized  compounds. 

•  Gene  expression  of  CCR5  and  CCL5  was  quantitated  by  SYBR-based  Real-time 
PCR  for  three  cell  lines,  P69,  Ml 2,  and  F6. 

•  Antiproliferation  assay  to  evaluate  the  anti  prostate  cancer  activity  of  the  new 
ligands  has  been  conducted  for  nine  new  ligands.  Three  different  prostate  cancer  cell 
lines  are  adopted  and  the  cells  are  treated  for  24,  48,  and  72  hours  with  varies 
concentrations  of  ligands  tested.  All  the  ligands  tested  showed  significant 
antiproliferation  activity. 

•  A  homology  modeling  study  has  been  conducted  to  characterize  the  3D 
conformation  of  chemokine  receptor  CCR5  and  its  interaction  with  the  lead  compound, 
anibamine. 
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Reportable  Outcomes 

One  manuscript  has  been  published  so  far  (see  Appendices)  and  two  more  manuscripts 
are  under  preparation  for  publication. 


Two  poster  presentations  related  to  this  project  have  been  presented  at  National  and 
International  level  symposiums.  They  are: 

1.  Joanna  L.  Adams,  Guo  Li,  Amanda  C.  Richardson,  Joy  L.  Ware,  and  Yan  Zhang. 
Prostate  Cancer  Cell  Line  Proliferation  Inhibited  by  Anibamine,  a  Natural  Product  CCR5 
Antagonist.  Summit  on  Systems  Biology  2007,  June  5-7,  2007,  Richmond,  Virginia. 

2.  Kendra  Haney,  Guo  Li,  John  Bajacan,  Yan  Zhang.  Natural  Product  Based  Novel  Anti- 
Prostate  Cancer  Drug  Discovery.  Massey  Cancer  Center  Research  Retreat,  Virginia 
Commonwealth  University,  June  2008,  Richmond,  Virginia. 


One  patent  application  has  been  prepared  and  the  provisional  invention  related  will  be 
fded  in  due  course. 


One  Master  Degree  student  (Joanna  Adams)  has  graduated  after  finishing  a  research 
project  related. 

One  Master  Degree  student  (Kendra  Haney)  is  expected  to  graduate  in  the  Summer  of 
2009  upon  the  completion  of  thesis  based  on  this  project. 


One  Postdoctoral  Associate  has  been  trained  under  the  support  of  this  award  and  is  ready 
to  apply  for  industrial  chemist  position. 
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Conclusion 


The  major  focus  of  the  research  pro  ject  will  be  the  syntheses  of  the  ligands  we  designed 
as  chem  okine  receptor  CCR5  antagonists  with  novel  structural  feature(s)  derived  from 
anibamine,  the  first  natural  product  as  a  CCR5  antagonist.  The  biological  characterization 
of  these  ligands  for  anti  prostate  cancer  activity  is  being  conducted.  Once  biological  test 
results  are  available,  Structure-Activity-Relationship  (SAR)  analyses  will  be  conducted  to 
identify  the  structural  feat  ures  for  the  next  generation  lead  compound(s).  When  enough 
data  is  available,  3D  QSAR  studies  will  be  performed.  Concomitantly,  we  will  refine  our 
already  completed  homology  model  for  the  CCR5  receptor  in  response  to  these  data. 

So  far  we  have  finished  the  chemical  synt  hesis  of  thirty-two  novel  ligands  that  are 
designed  based  on  “Deconstruction-Reconstr  uction-Elaboration”  method.  A  RANTES- 
binding  inhibition  assay  protocol  has  been  set  u  p  and  ready  to  test  our  drug  candidates. 
Other  biological  screening  includes  the  te  st  of  capacity  of  these  novel  com  pounds  to 
inhibit  proliferation  and/or  apoptosis  by  the  human  prostate  cancer  cell  lines  Ml 2,  PC-3, 
DU- 145  and  LNCaP  has  been  co  nducted  to  ev  aluate  th  eir  efficacy.  Also  we  have 
assessed  changes  in  the  expression  of  CCL5  and  CCR5  in  hu  man  prostate  epithelial  cell 
lines  of  the  sa  me  cellular  lineage  but  with  differing  in  vivo  phenotypes  (P69,  F6,  and 
M12).  A  molecular  m  odeling  study  has  been  conducted  to  elucidate  the  interaction  of 
anibamine,  the  lead  compound,  to  the  targ  et  protein,  chemokine  receptor  CCR5,  to 
further  facilitate  our  next  generation  molecular  design. 

In  the  nex  t  year,  we  will  f  ocus  on  the  stru  cture  m  odification/elaboration  pa  rt  of  the 
proposal  by  introducin  g  m  ore  varieties  of  substitutions  at  th  e  two  side  chain  s  by 
following  the  “Lipinski’s  rule  of  5”  fo  r  dr  ug-like  compounds.  Further  biological 
screening  will  start  from  radioligand  binding  assay,  to  in  vitro  an  ti  prostate  cancer  tests. 
We  wish  to  identify  a  prom  ising  lead  co  mpound  with  novel  chem  ical  and  structural 
features  as  well  as  high  therapeutic  indices  to  serve  our  long-term  goal  in  finding  a  potent 
anti  prostate  cancer  agent. 

The  medical  product  expected  from  this  project  in  long  term  is  discovery  of  a  novel  and 
potent  anti  prostate  cancer  agent  with  new  therapeutic  index. 
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ABSTRACT 


i=> 

1 

Anibamine 

The  total  synthesis  of  anibamine,  the  first  and  only  natural  product  known  as  a  chemokine  receptor  CCR5  antagonist,  is  reported  herein. 
Anibamine  was  synthesized  from  acetylacetone  and  cyanoacetamide  in  10  steps. 


.  0  0  o 


Anibamine,  a  novel  pyridine  quaternary  alkaloid  recently 
isolated  from  Aniba  sp.,12  has  been  found  to  effectively  bind 
to  the  chemokine  receptor  CCR5  with  an  IC50  at  1  /uM  in 
competition  with  125I-gpl20,  a  HIV  viral  envelop  protein 
binding  to  CCR5  with  high  affinity.1  The  chemokine  receptor 
CCR5,  a  G-protein-coupled  receptor,  has  been  identified  as 
an  essential  co-receptor  for  HIV  virus  entry  to  host  cells.3 
Therefore,  an  antagonist  of  CCR5  receptor  that  would  inhibit 
the  cellular  entry  of  human  immunodeficiency  virus  type  I 
(HIV-1)  provides  a  new  therapy  choice  for  the  treatment  of 
HIV  infection.4 

Currently,  all  the  known  antagonists  to  CCR5  have  been 
developed  by  optimization  of  the  lead  compounds  from  high- 
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throughput  screening,  such  as,  SCH-D,5  TAK220,6  and 
UK4278577  (Figure  1). 


UK-427857  (In  Clinical  Trial)  TAK  220  (Clinical  Candidate) 


Figure  1.  Anibamine  and  some  known  CCR5  antagonists 


Anibamine  is  the  first  and  the  only  natural  product  that 
has  been  identified  as  a  CCR5  antagonist  with  high  affinity. 
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T.  N.  C.  Nat.  Rev.  Drug  Discovery  2002, 1,  347—358.  (d)  Lusso,  P.  EMBO 
J.  2006,  25,  447-456. 
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Clearly,  anibamine  possesses  a  novel  structural  skeleton 
compared  with  all  other  known  CCR5  antagonists.5  Thus  the 
chemical  synthesis  of  anibamine  and  its  analogues  may  lead 
to  a  new  class  of  AIDS  therapeutic  agents.  Herein  is  the 
first  report  of  the  total  synthesis  of  anibamine. 

As  shown  below,  anibamine  has  a  unique  structural 
skeleton,  in  which  a  pyridine  quaternary  ion  is  embodied  by 
a  fused  five-member  ring.  Such  a  skeleton  can  also  be  called 
the  2,3-dihydro-7/7-indolizinium  ion.  On  the  pyridine  ring, 
two  ten-carbon  aliphatic  chains  are  attached  at  positions  3 
and  5  via  cis  double  bonds.  In  addition,  both  positions  2 
and  4  are  occupied  by  methyl  groups.  The  retrosynthetic 
analysis  of  anibamine  is  illustrated  in  Scheme  1. 


Scheme  1.  The  Retrosynthetic  Analysis  of  Anibamine 


The  two  ten-carbon  side  chains  can  be  introduced  by 
Wittig  reactions  between  the  Wittig  reagent  of  1-bromo- 
nonane  and  the  3, 5 -dialdehyde  intermediate  2.  Intermediate 

2  can  be  prepared  by  the  reduction  of  the  3,5-dicyano 
intermediate.  The  fused  five-member-ring  system  can  be 
achieved  by  introducing  a  three-carbon  side  chain  at  position 
2,  followed  by  the  cyclization  reaction  to  form  the  indoli- 
zinium  ring.  Apparently  the  substituted  pyridine  3  is  the  most 
critical  intermediate  since  it  possesses  all  the  potential 
functional  groups. 

Two  synthetic  routes  were  explored  to  prepare  the 
intermediate  3,  as  shown  in  Scheme  2.  In  the  first  route,  the 
condensation  reaction  between  acetylacetone  and  N,N- 
dimethylformide  dimethyl  acetal  provided  compound  4. 8 
Then  4  was  refluxed  with  hydroamine  hydrochloride  in 
methanol  to  give  compound  5. 9  Next  5  was  coupled  with 
aniline10  to  provide  compound  6.  Finally  the  cyclization 
reaction  between  6  and  malononitrile11  led  to  the  intermediate 

3  in  an  overall  30%  yield.  In  the  second  route,  using  the 


(5)  Palani,  A.;  Tagat,  J.  R.  J.  Med.  Chem.  2006,  49,  2851-2857. 

(6)  Imamura,  S.;  Ichikawa,  T.;  Nishikawa,  Y.;  Kanzaki,  N.;  Takashima, 
K.;  Niwa,  S.;  Iizawa,  Y.;  Baba,  M.;  Sugihara,  Y.  J.  Med.  Chem.  2006,  49, 
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Scheme  2.  The  Chemical  Synthesis  of  Intermediate  3 
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procedure  of  Haley  et  al.,12  the  condensation  reaction  between 
acetylacetone  and  cyanoacetamide  gave  compound  7,  which 
was  brominated  with  NBS  in  TFA  and  H2SO4  to  give  8.13 
Next,  the  nucleophilic  cyanization  of  compound  8  was 
achieved  by  refluxing  with  cuprous  cyanide  in  DMF  to  give 
3. 13  Here,  the  total  yield  was  46%  over  3  steps.  The  first 
route  is  one  step  longer  than  the  second  one,  yet  is  more 
convenient  for  large-scale  synthesis.  Since  two  cyano  groups 
can  be  introduced  subsequently  in  the  second  route,  it  may 
provide  the  possibility  of  attaching  two  different  side  chains 
on  positions  3  and  5.  This  will  be  critical  to  anibamine 
analogue  synthesis  in  the  future. 

The  bromination14  of  3  with  TBAB/P205  led  to  the 
2-bromo-substituted  intermediate  9  (Scheme  3).  The  pal- 


Scheme  3.  The  Synthesis  of  Intermediate  11 


ladium-catalyzed  alkynylation  reaction  of  compound  9  with 
protected  propargyl  alcohol  gave  compound  10. 15,16 


(12)  Haley,  C.  A.  C.;  Maitland,  M.  J.  Chem.  Soc.  1951,  3155-3174. 

(13)  Mai,  P.;  Lourderaj,  U.;  Parveen;  Venugopalan,  P.;  Moorthy,  N.  J.; 
Sathyamurthy,  N.  J.  Org.  Chem.  2003,  68,  3446-3453. 

(14)  Kato,  Y.;  Okada,  S.;  Tomimoto,  K.;  Mase,  T.  Tetrahedron  Lett. 
2001,  42,  4849-4851. 

(15)  Lautens,  M.;  Yoshida,  M.  J.  Org.  Chem.  2003,  68,  762-769. 

(16)  Negishi,  E.;  Anastasia,  L.  Chem.  Rev.  2003,  103,  1979-2017. 
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To  our  delight,  the  hydrogenation  of  compound  10  under 
Pd/C  in  methanol  at  room  temperature  gave  intermediate  11 
in  quantitative  yield.  After  that,  several  different  reducing 
agents  were  applied  to  prepare  intermediate  2.  These  results 
are  summarized  in  Table  1. 


Table  1.  Preparation  of  Intermediate  2 


11  2  2'  2" 


yield 


entry 

reduction  condition 

product 

(%) 

1 

Raney-Ni,  NaH2P02,  Py/AcOH/H20 

no  reaction 

2 

Raney-Ni,  NaH2P02,  AcOH/H20 

2' 

70 

3 

Raney-Ni,  5  bar  of  H2,  AcOH 

2' 

95 

4 

Raney-Ni,  1  bar  of  H2,  AcOH 

2 

90 

5 

Raney-Ni,  75%  HCOOH  aq 

2 " 

13 

6 

DIBAL-H,  THF,  0  °C  to  reflux 

no  reaction 

7 

DIBAL-H,  toluene,  0  °C 

2 

77 

With  use  of  Raney-Ni17’18  as  the  reduction  reagent  in 
different  solvent  systems,  the  3, 5 -diaminomethyl- substituted 
product  (2')  was  the  only  one  except  for  the  procedure  with 
HCOOH  as  solvent,  which  led  to  compound  2"  in  13%  yield. 
When  DIBAL-H  was  used  as  the  reduction  reagent  in 
toluene,  the  desired  product  (2)  was  obtained  in  77% 
separation  yield,19  in  contrast  to  the  unsuccessful  reaction 
in  THF.  Our  explanation  of  this  result  is  that  the  activity  of 
the  electrophilic  reduction  reagent  DIBAL-H  in  toluene  is 
higher  than  that  in  the  comparatively  more  polar  solvent 
THF.2021 

To  introduce  the  aliphatic  side  chains  onto  positions  3  and 
5,  the  Wittig  reaction  was  explored  under  different  basic 
conditions,  such  as  NaH/toluene,22  NaH/DMSO,23  LHMDS,24 
and  n- BuLi.25  No  reaction  was  observed  when  NaH  was  used 
as  the  base  in  toluene.  A  trace  amount  of  the  target  product 
was  observed  while  most  of  the  starting  material  decomposed 
under  NaH/DMSO  and  n-BuLi  conditions.  Finally,  interme¬ 
diate  12  was  obtained  as  the  major  product  when  LHMDS 
was  adopted  as  the  base.  We  noticed  that  the  crude 


(17)  Rapold,  T.;  Senn,  M.  U.S.  Patent,  5,384,403,  1995. 

(18)  Heffner,  R.  J.;  Jiang,  J.  J.;  JoulliG,  M.  M.  J.  Am.  Chem.  Soc.  1992, 
114,  10181-10189. 

(19)  Song,  Z.  J.;  Zhao,  M.;  Desmand,  R.;  Devine,  P.;  Tschaen,  D.  M.; 
Tillyer,  R.;  Frey,  L.;  Heid,  R.;  Xu,  F.;  Foster,  B.;  Li,  J.;  Reamer,  R.;  Volante, 
R.;  Grabowski,  E.  J.  J.;  Dolling,  U.  H.;  Reider,  P.  J.  J.  Org.  Chem.  1999, 
64,  9658-9667. 
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1448.  (b)  Heathcock,  C.  H.;  Finkelstein,  B.  L;  Jarvi,  E.  T.;  Radel,  P.  A.; 
Hadley,  C.  R.  J.  Org.  Chem.  1988,  53,  1922-1942. 


dialdehyde  2  can  be  used  in  this  step  without  further 
purification  after  the  DIBAL-H  reduction  reaction,  and  the 
combined  yield  was  71%  over  two  steps  from  compound 

11. 

A11  four  possible  isomers  were  observed  as  products  of 
this  reaction  (Scheme  4).  Various  chromatographic  conditions 


Scheme  4.  The  Wittig  Reaction 


38% 


have  been  applied  yet  no  successful  separation  of  any  one 
of  these  four  isomers  has  occurred.  Therefore,  after  primary 
purification,  the  mixture  was  used  for  the  next  step.  By 
calculating  the  ratios  of  the  integrates  of  the  methyl  group 
at  position  6  on  the  pyridine  ring  in  the  !H  NMR  spectrum, 
the  composition  of  these  four  isomers  in  the  mixture  were 
determined,  and  compound  12  was  characterized  as  the  major 
one. 

Interestingly,  the  most  common  PMB  deprotection  reagent, 
DDQ,  was  proved  inefficient  in  both  CH2CI2/H2O26’27  and 
CH2CI2/  buffer  (pH  7)  system28  for  compound  12.  In  fact, 
the  PMB  group  was  removed  quickly  with  reasonable  yield 
under  acidic  conditions  to  give  compound  13. 29 

The  five-member-ring  closure  was  achieved  by  treating 
compound  13  with  methanesulfonyl  chloride  and  triethyl- 
amine  at  room  temperature.30’31  The  crude  product  was 
purified  by  preparative  HPLC  with  the  previously  reported 
condition.1  Both  anibamine  (1)  and  its  (IIE,22E)  isomer  were 
isolated  as  the  trifluoroacetic  acid  salts.  The  spectral  proper¬ 
ties  of  anibamine  were  compared  with  those  in  the  literature 
and  no  significant  differences  were  observed. 


Scheme  5.  The  Synthesis  of  Anibamine 


The  anti-HIV  activity  of  anibamine  and  its  (IIE,22E) 
isomer  was  evaluated.  In  an  assay  that  determines  inhibition 
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of  HIV-lBaL  attachment  to  GHOST  R5  cells,  anibamine 
showed  an  EC50  at  0.6  juM  and  its  (l\E,22E)  isomer  did  at 
0.8  /uM.32 

In  summary,  anibamine  was  synthesized  from  acetyl- 
acetone  and  cyanoacetamide  in  10  steps  with  7.9%  overall 
yield,  or  from  acetylacetone  and  A,A-dimethylformide  di- 


(26)  Oikawa,  Y.;  Yoshioka,  T.;Yonemitsu,  O.  Tetrahedron  Lett.  1982, 
23,  885-888. 

(27)  Horita,  K.;  Yoshioka,  T.;  Tanaka,  T.;  Oikawa,  Y.;  Yonemitsu,  O. 
Tetrahedron  1986,  42,  3021—3028. 

(28)  Yu,  W.  S.;  Zhang,  Y.;  Jin,  Z.  D.  Org.  Lett.  2001,  3,  1447-1450. 

(29)  Jenkins,  D.  J.;  Riley,  A.  M.;  Potter,  B.  V.  L.  J.  Org.  Chem.  1996, 
61,  7719-7726. 

(30)  Ciufolini,  M.  A.;  Roschangar,  F.  J.  Am.  Chem.  Soc.  1996,  118, 
12082-12089. 

(31)  Bates,  R.  W.;  Boonsombat,  J.  J.  Chem.  Soc.,  Perkin  Trans.  2001, 
1,  654-656. 

(32)  Dextran  sulfate  were  provided  as  a  relevant  positive  control 
compound  for  the  individual  assays,  with  anti-HIV  activity  of  EC50  at  14.5 
fiM. 


methyl  acetal  in  11  steps  with  5.1%  overall  yield.  To  our 
knowledge,  this  is  the  first  report  of  the  total  synthesis  of 
anibamine.  The  above  synthetic  routes  also  offer  the  op¬ 
portunity  to  prepare  anibamine  analogues  for  further  bio¬ 
logical  evaluation  and  SAR  study  for  their  anti-HIV 
activity. 
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Anibamine,  a  novel  pyridine  quaternary  alkaloid  recently  isolated  from  Aniba  sp .,  has  been  found  to  effectively 
bind  to  the  chemokine  receptor  CCR5  with  an  IC50  at  1  pM  in  competition  with  125I-gpl20,  an  HIV  viral 
envelope  protein  binding  to  CCR5  with  high  affinity.  Since  CCR5,  a  G-protein-coupled  receptor,  is  an 
essential  coreceptor  for  the  human  immunodeficiency  virus  type  I  (HIV-1)  entry  to  host  cells,  a  CCR5 
antagonist  that  inhibits  the  cellular  entry  of  HIV-1  provides  a  new  therapy  choice  for  the  treatment  of  HIV. 
Anibamine  provides  a  novel  structural  skeleton  that  is  remarkably  different  from  all  lead  compounds  previously 
identified  as  CCR5  antagonists.  Here,  we  report  comparative  docking  studies  of  anibamine  with  several 
other  known  CCR5  antagonists  in  two  CCR5  homology  models  built  based  on  the  crystal  structures  of 
bovine  rhodopsin  and  human  /^-adrenergic  receptor.  The  binding  pocket  of  anibamine  has  some  common 
features  shared  with  other  high  affinity  CCR5  antagonists,  suggesting  that  they  may  bind  in  similar  binding 
sites  and/or  modes.  At  the  same  time,  several  unique  binding  features  of  anibamine  were  identified,  and  it 
will  likely  prove  beneficial  in  future  molecular  design  of  novel  CCR5  antagonists  based  on  the  anibamine 
scaffold. 


INTRODUCTION 

According  to  UNAIDS  Report  2006, 1  an  estimated  38.6 
million  people  worldwide  were  living  with  HIV  at  the  end 
of  2005.  In  the  United  States  about  one  million  people  are 
currently  living  with  HIV.2  Since  the  AIDS  pandemic 
remains  one  of  the  leading  global  public  health  threats,  the 
continued  discovery  and  development  of  new  antiretroviral 
drugs  with  reduced  toxicity,  enhanced  potency,  novel  mech¬ 
anisms  of  action,  and  reduced  prevalence  of  adverse  drug- 
drug  interactions  remains  a  very  high  priority.  The  chemokine 
receptors  CCR5  and  CXCR4  are  essential  coreceptors 
required  for  the  attachment  of  HIV-1  to  the  host  cell 
membrane3-5  suggesting  that  CCR5  antagonists  may  lead 
to  a  novel  therapeutic  strategy  in  blocking  the  entry  of  HIV 
into  host  cells,  especially  when  utilized  in  combination  with 
the  classical  reverse  transcriptase  and  protease  inhibitors.6,7 
At  the  same  time,  the  CCR5-A32  homo-  and  heterozygous 
genotype  carriers  have  been  shown  to  confer  either  resistance 
to  HIV-1  infection  or  delayed  progression  of  the  disease.8,9 
Therefore,  it  is  clear  that  CCR5  plays  an  essential  role  in 
HIV  pathogenesis. 

The  chemokine  receptor  CCR5  (Figure  1)  belongs  to  the 
extensively  exploited  family  of  G-protein-coupled  receptors 
(GPCRs).  Chemokines  and  chemokine  receptors  play  a 
crucial  role  in  the  trafficking  of  leukocyte  populations  in  the 
body  and  are  involved  in  the  development  of  a  wide  variety 
of  human  diseases.10  Despite  the  fact  that  it  has  been  more 
than  ten  years  since  the  chemokine  receptor  CCR5  was 
identified  as  one  of  the  essential  coreceptors  utilized  by 
HIV-1  for  entry  into  human  CD4+  cells,3-5  there  are 
currently  only  a  few  CCR5  antagonists  (Figure  2)  being 
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investigated  as  anti-HIV- 1  agents  in  human  clinical  trials.11-14 
In  fact,  only  one  drug  (Maraviroc  or  Selzentry)  was  approved 
by  the  FDA  on  August  6,  2007 14d  even  though  the  FDA 
raised  concerns  that  Maraviroc  could  be  associated  with  an 
increased  risk  of  liver  damage,  lymphoma,  infections,  and 
increased  risk  of  heart  attack. 14c,d  Also,  clinical  trials  of 
Aplaviroc  were  halted  due  to  severe  liver  side  effects  in  both 
treatment-naive  and  treatment-experienced  patients.136  It 
should  be  noted  that  the  lead  compounds  for  all  of  these 
CCR5  antagonists  were  discovered  through  high-throughput 
screening  campaigns  targeting  the  disruption  of  the  virus 
attachment  to  the  coreceptor.  While  this  is  often  a  valid 
approach  that  yields  interesting  lead  compounds,  natural 
products  can  often  be  more  diverse,  having  novel  structure 
skeletons  with  wider  ranges  of  chemical  shape  and  structural 
features.  It  is  well-known15  that  natural  products  offer 
unmatched  structural  variety  compared  with  combinatorial 
chemistry  libraries,  and  their  high  diversity  and  specific 
biological  activities  very  often  make  them  useful  resources 
for  drug  discovery. 

Recently,  anibamine,  a  novel  pyridine  quaternary  alkaloid 
that  was  isolated  from  Aniba  sp .,  has  been  found  to 
effectively  bind  to  CCR5  with  an  IC50  of  1  juM  in  competitive 
binding  assays  with  125I-gpl20,  a  viral  envelope  protein  that 
binds  to  CCR5  with  high  affinity.16  As  the  first  and  only 
natural  product  to  date  showing  CCR5  antagonism,  anibam¬ 
ine  has  a  novel  structure  quite  different  from  those  other 
known  classes  of  CCR5  antagonists.  Recently,  we  reported 
the  total  synthesis  of  anibamine17  applying  a  flexible 
synthetic  scheme  that  will  allow  the  production  of  a  wide 
variety  of  analogs. 

It  should  be  noted  that  the  1  juM  activity  of  anibamine  is 
in  line  with  those  of  other  “lead”  compounds  (0.39—2  pM) 
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Figure  2.  Chemical  structures  of  anibamine  and  several  known  potential  CCR5  antagonists. 


that  were  developed  into  the  antagonists  illustrated  in  Figure 
2. 17-21  Anibamine  also  has  structural  similarities  to  these 
compounds,  in  particular  a  potentially  charged  nitrogenic 
central  moiety  with  two  or  more  hydrophobic  side  chains 
attached.  However,  because  of  its  structural  novelty,  we 
cannot  dismiss  the  possibility  that  anibamine  may  possess  a 
different  binding  modality  than  that  proposed  for  the  other 
CCR5  antagonists.22  In  the  absence  of  experimental  structural 
data  for  the  CCR5  receptor,  one  of  the  few  ways  to 
investigate  the  binding  modes  of  these  compounds  is  through 
molecular  modeling.  An  understanding  of  anibamine  and  its 
binding  mode  in  CCR5  may  assist  the  design  of  novel 
derivatives  that  are  CCR5  antagonists  with  higher  affinity. 
Since  anibamine  has  a  unique  chemical  structural  skeleton 
with  a  potentially  novel  molecular  design  template,  compar¬ 
ing  its  binding  modality  with  those  of  other  known  CCR5 
antagonists  should  be  instructive  regarding  its  activity  and 
further  validate  it  as  a  lead  compound  for  next  generation 
molecular  design. 

In  general,  computational  and  mathematical  approaches, 
such  as  structural  bioinformatics,23  molecular  docking,24 
molecular  packing,25  pharmacophore  modeling,26  Monte 


Carlo  simulated  annealing  approach,27  graph/diagram  ap¬ 
proach,28  diffusion-controlled  reaction  simulation,29  biom- 
acromolecular  internal  collective  motion  simulation,30 
QSAR,31  protein  subcellular  location  prediction,32  identifica¬ 
tion  of  membrane  proteins  and  their  types,33  identification 
of  enzymes  and  their  functional  classes,34  identification  of 
GPCR  and  their  types,35  identification  of  proteases  and  their 
types,36  protein  cleavage  site  prediction,37  and  signal  peptide 
prediction38  can  provide  useful  information  and  insights  for 
both  basic  research  and  drug  design  and  hence  are  widely 
welcome  by  the  science  community. 

More  specifically,  in  the  absence  of  specific  structural  data, 
molecular  modeling  studies  that  build  and  refine  homology 
models  of  the  receptor  and  dock  the  various  ligands  into  these 
models  are  the  most  accessible  approaches  to  studying 
binding  modes  for  ligands.  Like  nearly  all  GPCRs,  a  crystal 
structure  for  the  CCR5  receptor  is  not  currently  available, 
and  a  template-derived  homology  model  has  to  be  used  to 
gain  insight  into  the  receptor  structure.  In  fact,  a  homology 
modeling  method  based  on  the  crystal  structure  of  bovine 
rhodopsin39  has  been  successfully  applied  to  other  GPCRs,40 
including  the  CCR5  receptor,22  to  further  understand  the 
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ligand-protein  interactions  and  to  identify  new  and  potent 
ligands.  For  example,  a  homology  model  of  Angiotensin  II 
type  1  (ATI)  receptor  was  used  to  explore  the  binding  sites 
of  several  nonpeptide  ATI  receptor  antagonists,41  and  a 
homology  model  of  Ml  muscarinic  acetylcholine  receptor 
was  applied  to  understand  the  mechanism  by  which  the 
agonist-receptor  complex  activates  the  G  proteins.42  More 
relevant,  the  models  of  CCR5  developed  by  several  groups 
helped  characterize  the  binding  mode  of  known  CCR5 
agonists  and  antagonists.22  It  is  believed  that  with  all  the 
lessons  learned  from  previous  experience,  homology  model¬ 
ing  of  GPCRs  based  on  bovine  rhodopsin  crystal  structure 
will  help  in  structure-based  drug  design  and  virtual  screening 
for  therapeutic  applications  at  many  GPCR  targets.  The 
recent  seminal  resolution  of  two  human  /?2-adrenergic 
receptor  (/32AR)  crystal  structures43  led  us  to  question  which 
crystal  structure  would  be  the  more  reasonable  template  for 
building  the  homology  model  of  the  CCR5  receptor.  Both 
/?2AR  structures  have  a  bound  inverse  agonist  carazolol,  and 
presumably  these  structures  were  closer  to  the  inactive  state 
of  /32AR.  However,  the  broken  ionic  locks  in  the  reported 
structures  and  the  increased  affinity  of  /?2 AR-T4  chimera  for 
agonists  compared  to  the  wild  type  /32AR  suggest  that  these 
/32AR  crystal  structures  seemed  closer  to  the  active  state  of 
the  /?2AR  receptor  or,  possibly,  a  state  between  active  and 
inactive.44  Thus,  a  homology  model  based  on  the  /J2AR 
structure  would  also  be  worth  exploring  for  valuable 
information  of  CCR5  receptor. 

In  this  study,  the  binding  mode  of  anibamine  in  the 
antagonist  locus  of  CCR5  has  been  studied  by  building  and 
analyzing  two  complete  (i.e.,  including  both  the  extracellular 
and  intracellular  loops)  homology  models  of  the  CCR5 
receptor.  As  in  nearly  all  homology  models  of  GPCRs,  one 
of  our  models  was  based  on  the  crystal  structure  of  bovine 
rhodopsin  in  the  dark  state  (PDB  code  1F88),  which  was, 
up  until  very  recently,  the  only  available  crystal  structure 
for  a  GPCR.  The  other  homology  model  was  constructed 
based  on  the  crystal  structure  of  /32AR  (PDB  code  2RH1). 
The  docking  studies  of  known  CCR5  antagonists  into  each 
of  these  two  CCR5  homology  models  showed  that  the  CCR5 
homology  model  based  on  the  bovine  rhodopsin  crystal 
structure  seemed  more  compatible  with  the  site-directed 
mutagenesis  study  results.  The  docking  studies  of  anibamine 
in  these  two  homology  models  showed  that  the  anibamine 
binding  mode  shares  some  common  binding  features  with 
other  CCR5  antagonists.  Interestingly,  some  unique  inter- 
molecular  interactions  between  anibamine  and  the  CCR5 
receptor  homology  model  may  further  direct  structural 
modifications  of  anibamine  as  a  lead  for  next  generation 
CCR5  antagonist  design. 

RESULTS  AND  DISCUSSION 

The  computational  analysis  of  anibamine  was  processed 
in  several  steps:  first,  we  performed  a  conformational  study 
of  it  and  other  known  CCR5  antagonists.  Next,  homology 
models  of  the  chemokine  receptor  CCR5,  including  trans¬ 
membrane,  extracellular  and  intracellular  components,  were 
built  based  on  the  crystal  structure  of  bovine  rhodopsin  and 
human  /32AR,  followed  by  molecular  dynamics  simulation/ 
annealing  to  optimize  the  conformations  of  the  models. 
Finally,  anibamine  and  the  other  ligands  were  docked  into 


these  two  CCR5  models,  with  particular  emphasis  on 
evaluating  the  CCR5  binding  pocket  interactions  with 
anibamine  and  other  known  CCR5  antagonists. 

Conformation  Study  of  Anibamine  and  Other  Known 
CCR5  Antagonists.  In  order  to  simulate  the  conformation 
of  anibamine  in  a  biological  system  as  well  as  in  the  CCR5 
binding  site,  a  conformation  study  was  conducted  by 
minimization  and  molecular  dynamics  simulations.  Other 
compounds  reported  to  be  potential  CCR5  antagonists, 
including  Aplaviroc,  Maraviroc,  and  Vicriviroc  that  have 
entered  clinical  trials,  were  similarly  analyzed.  The  lowest 
energy  conformation  obtained  for  each  ligand  was  further 
simulated  in  a  solvated  system,  and  these  results  were  then 
applied  to  the  initial  conformation  to  explore  their  binding 
modes  in  the  CCR5  receptor  antagonist  binding  locus.  These 
conformations  were  only  expected  to  provide  a  reasonable 
starting  point  for  our  docking  studies;  they  were  to  endure 
some  possibly  dramatic  alterations  during  docking  in  order 
to  achieve  the  lowest  energy  conformations  for  the  ligand— receptor 
complex  (vide  infra). 

Homology  Modeling  of  Chemokine  Receptor  CCR5. 

Chemokine  receptors  belong  to  the  rhodopsin  subfamily  in 
the  GPCR  superfamily  and  are  characterized  by  a  heptapep- 
tidic  transmembrane  helical  fold  (7-TM)  spanning  the  plasma 
membrane  (see  Figure  1).  Homology  modeling  studies  of 
the  CCR5  receptor,  based  on  the  crystal  structures  of  the 
dark  state  bovine  rhodopsin  at  2.8  A  resolution  (PDB  code 
lF88)39b  and  on  the  human  /32AR  structure  at  2.4  A 
resolution  (PDB  code  2RHl),43a  involved  several  steps: 
sequence  alignment,  assignment  of  Cartesian  coordinates  for 
amino  acid  residues  in  structurally  conserved  regions, 
assignment  of  coordinates  for  amino  acid  residues  in 
structurally  nonconserved  regions,  and  refinement  of  coor¬ 
dinates  of  side  chains  of  the  entire  protein. 

The  initial,  and  likely  most  essential,  step  in  homology 
modeling  is  sequence  alignment  among  the  CCR5,  the  bovine 
rhodopsin,  and  the  human  /32AR.  This  was  performed  for 
the  entire  protein  sequences  by  using  the  Insightll/Homology 
module  (Figure  3).  In  order  to  confirm  the  alignment 
automatically  generated  by  the  software,  manual  alignment 
of  the  amino  acid  residues  was  also  performed  by  comparing 
the  chemical  structures  of  the  autoaligned  residues.  While 
the  entire  protein  sequences  homology  level  was  ap¬ 
proximately  20%  among  all  three  proteins,  their  transmem¬ 
brane  domains  gave  58%  homology  identity  between  the 
bovine  rhodopsin  and  CCR5  while  56%  between  human 
/?2AR  and  CCR5  (Table  1).  This  further  validated  applying 
these  two  crystal  structures  as  templates  to  construct  homol¬ 
ogy  models  of  the  CCR5  receptor.  It  is  important  to  note 
that  CCR5  receptor’s  sequence  contains  all  of  the  conserved 
amino  acid  residues  found  in  most  GPCRs,  and  this  provides 
an  important  guideline  for  sequence  alignment. 

The  overlapping  amino  acid  sequences  between  rhodopsin, 
/?2AR,  and  the  CCR5  receptor,  i.e.,  structurally  conserved 
regions,  were  divided  into  different  lengths  separated  by  gaps 
(most  of  which  are  in  the  loop  domains).  The  coordinates 
for  the  amino  acid  residues  in  these  regions  of  rhodopsin  or 
/32AR  were  assigned  to  the  corresponding  residues  in  the 
CCR5  receptor. 

The  remaining  amino  acid  residues  of  the  CCR5  receptor 
constitute  the  structurally  nonconserved  regions,  while  most 
of  them  are  in  the  extracellular  and  intracellular  loop  regions. 
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Figure  3.  Sequence  alignment  of  human  CCR5  (CCR5),  bovine  rhodopsin  (BOVN)  (PDB  code  1F88),  and  human  /?2AR  (B2AR)  (PDB 
code  2RH1).  The  Ballesteros-Weinstein  numbering  system  was  adopted  to  mark  all  the  conserved  amino  acid  residues  among  most  of  the 
GPCRs  and  colored  in  red.  The  CCR5  protein  was  numbered  accordingly  under  its  sequence.  The  secondary  structure  of  the  CCR5  receptor 
3D  conformation  based  on  bovine  rhodopsin  crystal  structure  was  marked  out  below  the  sequence. 


Table  1.  Sequence  Homology  Analysis  of  the  CCR5  Receptor, 
Bovine  Rhodopsin,  and  the  /^-Adrenergic  Receptor 


CCR5 

homology  level  % 

TM1 

TM2 

TM3 

TM4 

TM5 

TM6 

TM7 

bovine  rhodopsin 
(1F88) 

53 

63 

61 

57 

65 

58 

43 

/?2-adrenergic 
receptor  (2RH1) 

50 

50 

61 

52 

65 

58 

57 

The  coordinates  for  residues  in  these  regions  were  obtained 
by  searching  for  these  sequences  in  the  PDB.  The  adopted 
conformations  of  these  nonconserved  regions  all  have  close 
similarity  to  the  backbone  orientation  of  the  corresponding 
regions  in  rhodopsin  or  /^AR. 

The  side  chain  conformational  energies  for  the  entire 
protein  were  minimized,  while  the  backbone  was  constrained 


followed  by  a  vigorous  dynamics  simulation  to  reach  an 
energetically  reasonable  conformation.  Then  the  local  ge¬ 
ometry  of  energy-minimized  structure  was  checked  using  the 
Insightll/Homology  program.  All  amino  acid  residues  have 
reasonable  bond  lengths  and  bond  angles.  The  analysis  of 
cp,  Xu  and  angles  of  the  resulting  protein  conformations 
was  further  conducted  with  Procheck  4.1, 45  and  the  results 
were  shown  in  Figure  4  and  Table  2.  While  the  model 
generated  based  on  the  /^AR  template  has  better  parameters 
that  that  using  bovine  rhodopsin  template,  it  mainly  is 
because  of  the  parameter  difference  from  the  two  templates, 
rather  than  the  modeling  operation  (see  the  Supporting 
Information). 

The  most  significant  difference  between  the  crystal 
structures  of  bovine  rhodopsin  and  /^AR  is  the  short 
a-helical  structure  of  /32AR  in  the  EL2  loop  compared  to 
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Figure  4.  Ramachandran  plots  of  A)  rhodopsin-based  and  B)  /32 AR- 
based  CCR5  receptor  homology  models. 

Table  2.  Data  from  the  Ramachandran  Plot  of  the  Rhodopsin-  and 
/32~ Adrenergic  Receptor-Based  CCR5  Homology  Models 

MFRa  AARh  GARC 

-  -  -  number  of 

model  n  %  n  %  n  %  subtotal^  Gly  and  Pro  totale 

Rho-based  184  69.7  74  28.0  6  2.3  264  24  288 

/32AR-based  221  83.7  38  14.4  5  1.9  264  24  288 

a  The  most  favored  regions  (A,  B,  and  L).  b  The  additional 

allowed  regions  (a,  b,  1,  and  p).  cThe  generously  allowed  regions 
(-a,  -b,  -1,  and  -p).  rfThe  total  number  of  residues  evaluated  except 
glycines  and  prolines.  e  The  total  number  of  amino  acids  in  the 
receptor  models. 


the  /3-sheet  structure  in  bovine  rhodopsin.  For  most  GPCRs, 
the  EL2  loop  part  is  very  important  for  ligand  recognition 
with  their  receptor.46-54  Thus,  we  were  concerned  if  an 
a-helix  structure  in  the  EL2  of  CCR5  receptor  is  a  reasonable 
conformation.  This  possibility  (forming  an  a-helix  in  the  EL2 
loop  of  the  CCR5  receptor)  was  checked  with  nnpredict 
(available  online  from  the  Cohen  group  at  Department  of 
Cellular  and  Molecular  Pharmacology,  University  of  Cali¬ 
fornia,  San  Francisco)  by  predicting  the  secondary  structure 


Sequence  IS2AR_E2  MHWYRATHQEA  I  NCYANETCCDFFT 
Prediction  - EE - H  H  H  H  H  H  H - 

Sequence  CCR5_E2  FTRSQKEGLHYTCSSHF  PYSQYQFW 
Prediction  - H  H - E  EEE - 

Figure  5.  /32AR  and  CCR5  EL2  secondary  structure  prediction  (H 
=  helix,  E  =  strand,  -  =  no  prediction). 


Figure  6.  Homology  models  of  CCR5  overlapped  in  ribbon.  The 
model  based  on  bovine  rhodopsin  crystal  structure  is  in  red,  while 
the  one  based  on  /32AR  is  in  green. 


of  the  EL2  loop.  Comparing  the  nnpredict  results  for  the 
/32AR  EL2  loop  with  that  of  the  CCR5  EL2  loop  suggests  a 
much  lower  probability  of  forming  an  a-helix  structure  in 
CCR5  (Figure  5).  Therefore,  this  domain  in  the  homology 
model  of  CCR5  constructed  from  the  /32AR  template  may 
need  to  be  further  scrutinized,  e.g.,  by  site-directed  mutagen¬ 
esis  data. 

The  differences  between  the  two  models  were  assessed 
(Figure  6).  The  overall  weighted  Root  Mean  Square  Distance 
(rmsd)  value  between  them  was  5.79  A.  In  contrast  to  the 
rhodopsin-derived  model,  the  TM1  domain  of  the  /32AR- 
derived  model  has  a  more  straight  extension.  The  TM2  to 
TM7  domains  of  the  two  models  basically  adopted  fairly 
similar  orientations.  As  expected  from  the  discussion  above, 
the  EL2  domains  were  significantly  different;  i.e.,  the 
rhodopsin-derived  model  presented  a  /3-sheet,  while  the 
/32AR-derived  model  presented  an  a-helix  conformation.  In 
addition  to  these  differences  that  largely  resulted  from  their 
parent  templates,  different  orientations  of  amino  acid  residues 
in  these  two  models  transmembrane  regions  were  also 
observed.  These  differences  may  result  in  each  specific 
binding  pocket  having  different  spatial  distributions  of 
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Table  3.  Top  GoldScores  for  Anibamine  and  Three  Other  Known 
Antagonists  Docked  in  Both  Homology  Models  of  CCR5 


fitness 

anibamine 

Aplaviroc 

Maraviroc 

Vicriviroc 

model  based  on 
bovine  rhodopsin 
(1F88) 

51.57 

66.71 

54.01 

48.83 

model  based  on 
@2-  adrenergic 
receptor  (2RH1) 

50.53 

47.64 

49.22 

43.61 

hydrophilic  and  hydrophobic  regions.  This  will  be  discussed 
with  the  docking  study  results  below. 

Binding  Site  Study  of  Anibamine  As  an  Antagonist 
to  CCR5.  The  lowest  energy  conformation  of  anibamine  and 
the  other  antagonists  (Aplaviroc,  Vicriviroc,  and  Maraviroc) 
from  the  conformation  study  were  first  docked  interactively 
into  the  antagonist  binding  locus  of  CCR5  at  the  upper  part 
of  the  transmembrane  region  of  the  receptor  followed  by 
minimization  and  dynamics  simulation  of  the  resulting 
CCR5-ligand  complex.  This  operation  was  a  necessity 
because  the  preliminary  automated  docking  studies  revealed 
no  reasonable  binding  loci  in  either  model  due  to  the 
stringently  steric  arrangement  of  the  amino  acid  residues  in 
the  potential  binding  pocket.  By  conducting  the  interactive 
docking,  the  ligand  binding  locus  in  the  receptor  was  thus 
identified  for  the  following  automated  docking  study.  The 
orientation  of  the  molecule  skeleton  in  the  binding  locus  was 
directed  by  the  following:  first,  the  putative  ionic  interaction 
between  the  tertiary  nitrogenic  group  in  the  ligand  and  the 
carboxylic  group  of  Glu283  on  TM  helix  7  in  the  receptor; 
and  second,  facing  the  hydrophobic  side  chain  portion  of 
the  ligand  toward  the  hydrophobic  TM  helices.  The  dockings 
were  validated  with  a  short  minimization  (500  iterations)  and 
dynamics  simulation  (1000  steps,  1000  fs  total)  followed  by 
a  more  rigorous  optimization  (minimization  of  5000  iterations 
and  a  dynamics  simulation  of  100,000  steps,  100  ps  total). 
In  both  processes,  the  backbone  of  the  receptor  was  fixed  to 
prevent  the  disruption  of  the  receptor  a-helical  bundle.  The 
lowest  energy  conformation  from  the  latter  simulation  was 
extracted  and  saved  for  analysis. 

The  automated  docking  study  was  conducted  by  using 
GOLD  3.1 55  (Genetic  Optimization  for  Ligand  Docking)  with 
most  default  parameters.  GOLD  permits  full  ligand  flexibility 
and  partial  protein  flexibility  (side  chain  and  backbone 
flexibility  for  up  to  ten  user-defined  residues).  The  protein 
model  used  for  docking  each  ligand  was  the  lowest  energy 
conformation  resulting  from  the  interactive  docking  studies. 
The  oxygen  atom  OE1  in  Glu283  defined  a  10  A  active  site 
radius  in  the  CCR5  receptor.  A  number  of  trials  were 
conducted,  and  the  GOLDScores  from  each  docking  trial 
were  recorded.  The  top  ten  solutions  from  the  final  trial  were 
saved.  The  ligand— receptor  complex  from  the  docking  result 
with  the  highest  GOLDScore  was  extracted  (Table  3).  It 
seemed  that  there  was  no  difference  in  the  binding  of 
anibamine  in  both  models  from  the  viewpoint  of  GOLD- 
S cores,  while  all  other  three  ligands  tended  to  bind  better 
with  relatively  higher  GOLDScores  in  the  homology  model 
of  CCR5  based  on  bovine  rhodopsin  than  the  one  based  on 

yS2AR. 

The  final  conformations  (from  the  GOLD  docking)  of  the 
receptor— ligand  complexes  are  illustrated  in  Figure  7,  and 
the  major  amino  acid  residues  in  the  binding  sites  of  the 
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eight  complexes  are  summarized  in  Table  4.  Analyses  of 
the  binding  pockets  of  anibamine  and  the  other  antagonists 
revealed  that  all  the  antagonists  bound  to  similar  hydrophobic 
pockets  in  the  CCR5  models,  while  they  may  occupy 
different  subloci  leading  to  the  differences  in  GOLDScores. 
For  all  four  ligands  studies,  the  binding  pockets  in  both 
models  were  built  mainly  by  TM2,  TM3,  TM5,  TM6,  and 
TM7.  Note  that  residue  Glu283  on  TM7  is  shared  by  all 
four  antagonists  in  both  models,  which  is  in  concordance 
with  the  existence  and  importance  of  a  potentially  positively 
charged  core  in  the  ligands. 

As  an  amphiphilic  molecule,  anibamine  carries  a  positive 
charge  and  two  highly  hydrophobic  side  chains.  Thus,  it  is 
not  surprising  that  its  binding  pocket  is  composed  mainly 
of  aliphatic  hydrophobic  amino  acid  residues  in  both  models, 
while  the  positively  charged  nitrogen  atom  is  in  the  vicinity 
of  Glu283  to  form  a  putative  salt  bridge.  The  GOLDScores 
for  each  model  are  not  significantly  different  (Table  3).  It 
should  be  noted  that  we  do  not  observe  significant  interac¬ 
tions  between  anibamine  and  the  extracellular  loop  II  (EL2) 
of  the  protein  although  Serl79  and  Seri 80  from  EL2  were 
in  the  vicinity  of  the  binding  locus  of  anibamine  in  both 
homology  models  of  CCR5.  Since  this  loop  is  believed  to 
be  involved  in  the  binding  of  gpl20  to  the  receptor,46-51 
this  observation  may  explain  why  anibamine’ s  inhibition 
affinity  to  CCR5  was  not  particularly  outstanding.  While  in 
the  rhodopsin-derived  model  (Figure  7A,C,E,G)  Tyrl08 
seemed  to  play  important  roles  in  the  binding  of  all  four 
ligands,  for  anibamine  it  probably  only  contributed  to  the 
hydrophobic  interactions  between  the  ligand  and  the  receptor. 
Residue  Tyr251  seemed  to  be  important  to  the  binding  of 
both  anibamine  and  Aplaviroc.  Some  amino  acid  residues 
were  unique  for  only  anibamine;  e.g.,  His289  not  only 
interacted  with  one  of  the  long  aliphatic  chains  in  the  ligand 
molecule  but  also  might  confer  cation-pi  and/or  pi-pi  stacking 
interaction  with  the  positively  charged  core  ring  system  (as 
His289  was  observed  in  the  vicinity).  In  the  /32AR-derived 
model  (Figure  7B),  the  two  long  aliphatic  side  chains  of 
anibamine  adopted  a  completely  different  orientation  than 
they  do  in  the  rhodopsin-derived  model,  which  lead  to  a 
different  set  of  amino  acid  residues  in  the  binding  pocket. 
Thus,  instead  of  Phell2,  Phel09  and  lie  164  were  part  of 
the  hydrophobic  pocket.  Also,  Trpl90  and  Met279  appeared 
to  play  some  role  in  the  recognition  of  anibamine  by  the 
CCR5  antagonist  binding  loci,  while  His289  was  not 
observed  in  the  binding  pocket.  Apparently  the  significant 
difference  from  the  orientation  of  the  two  aliphatic  side 
chains  in  the  molecule  provided  useful  information  for  future 
validation  of  the  binding  modes  as  well  as  for  the  anibamine 
derivative  design. 

In  agreement  with  a  very  comprehensive  site-directed 
mutagenesis  study,226  Aplaviroc  showed  very  strong  interac¬ 
tions  with  ECL2  in  the  rhodopsin-derived  model  (Figure  1C). 
For  example,  Seri  80  on  EL2  may  form  hydrogen  bond 
interaction  with  the  molecule.  In  agreement  with  the  same 
report,  residues  Tyrl08,  Phell2,  Tyr251,  and  Glu283  were 
all  observed  in  the  binding  pocket  of  Aplaviroc,  while 
Asn252  seemed  to  contribute  uniquely  to  the  binding  of  the 
molecule.  Besides,  an  ionic  or  hydrogen  bonding  interaction 
between  the  carboxylic  acid  moiety  of  the  molecule  and 
Lysl97  was  also  observed  in  this  model.  On  the  other  hand, 
we  did  not  observe  a  similar  interaction  between  this  acidic 
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Figure  7.  Binding  modes  of  the  CCR5  antagonists  in  the  homology  model  of  CCR5  constructed  based  on  bovine  rhodopsin  (red)  and 
human  /?2AR  (green).  The  amino  acid  residues  of  the  binding  pockets  were  as  shown  in  capped  sticks  and  the  ligands  were  in  ball  and  stick, 
while  the  receptor  models  were  in  ribbon.  ELs  and  part  of  TMs  were  omitted  for  clarity.  Oxygen  atoms  were  in  red,  fluorine  in  green,  and 
nitrogen  in  deep  blue.  Carbon  atoms  and  the  bonds  were  colored  according  to  the  ligands.  A  and  B,  anibamine  (orange);  C  and  D,  Aplaviroc 
(yellow);  E  and  F,  Maraviroc  (purple);  and  G  and  H,  Vicriviroc  (cyan). 

moiety  and  Lysl91  that  was  proposed  in  the  literature,  while  In  the  /^AR-derived  CCR5  model  (Figure  7D),  the  binding 

Lysl97  seemed  to  be  less  important  in  the  same  report.226  pocket  for  Aplaviroc  was  formed  by  a  totally  different  set 
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Table  4.  Major  Amino  Acid  Residues  in  the  Binding  Pockets  of  CCR5  Antagonists  in  the  Two  Homology  Models  of  the  Receptor" 

model  based  on  bovine  rhodopsin  (1F88)  model  based  on  /32~ adrenergic  receptor  (2RH1) 

anibamine  E283,  Y251,  H289,  1198,  Y108,  F112  E283,  Y251,  M279,  W190,  Y108,  F109,  1164 

Aplaviroc  E283,  Y251,  S180,  K197,  Y108 ,  F112,  W248,  N252  E283,  Y251,  Y108 ,  Q280,  N258,  L255,  W190 

Maraviroc  E283,  Y251,  1198,  Y108,  M279,  N252,  V199,  W248,  M287  E283,  Y251,  Y108,  M279,  Q280,  W190,  FI  12,  M287 

Vicriviroc  E283,  Y251,  Y108,  1198,  T284  E283,  Y251,  Y108,  L255,  T284,  Q280,  W190,  F109 

a  The  residues  in  Bold  are  consistent  with  the  site-directed  mutagenesis  report.  The  residues  in  Italics  are  in  contrary  to  the  site-directed 
mutagenesis  report.  The  residues  in  Bold/Italics  are  the  ones  that  were  reported  controversially  in  the  site-directed  mutagenesis  reports. 


of  amino  acid  residues,  including  Gln280,  Asn258,  Leu255, 
Trpl90,  Tyr251,  and  Tyrl08,  with  Glu283  as  the  only 
exception.  Among  these  residues,  Gln280,  Asn258,  and 
Leu255  have  not  been  tested  in  any  site-directed  mutagenesis 
studies,  while  Trpl90,  Tyr251,  and  Tyrl08  were  deemed  not 
essential  for  ligand  binding  to  CCR5  in  another  report.22f 
Besides,  Seri 80  was  not  observed  in  the  binding  locus  since 
EL2  was  not  involved  in  the  binding  pocket  formation  in 
this  model.  This  may  explain  the  lower  GOLDScore  of 
Aplaviroc  in  the  /32AR-derived  CCR5  model,  while  plausible 
hydrogen  bonding  between  Seri 80  and  the  Aplaviroc 
molecule  may  contribute  significantly  in  the  rhodopsin- 
derived  model. 

Maraviroc  showed  weak  interactions  with  the  EL2  region 
in  the  rhodopsin-derived  model  (Figure  7E).  Similar  to 
Vicriviroc,  Glu283,  Tyrl08,  and  lie  198  were  important  in 
the  binding  pocket  of  Maraviroc  in  this  model,  and  this  was 
also  supported  by  site-directed  mutagenesis  studies  22f  Also, 
Tyr251  was  part  of  the  binding  locus  in  this  model,  again  as 
supported  by  site-directed  mutagenesis. 22f  Moreover,  Trp248 
was  in  the  vicinity  of  the  ligand  molecule  in  this  model,  while 
it  was  not  supported  by  the  same  literature.  In  the  /32AR- 
derived  CCR5  model  (Figure  7F),  the  binding  pocket  for 
Maraviroc  showed  some  similarity  to  the  model  based  on 
rhodopsin,  e.g.  the  presence  of  Tyr251  and  Tyrl08  in  the 
binding  locus.  In  addition,  the  involvement  of  Phell2  and 
Trpl90  has  not  been  tested,  while  the  involvement  of  Met287 
in  the  binding  locus  was  not  supported  by  site-directed 
mutagenesis. 22f  While  the  Goldscores  in  the  two  models  did 
not  show  significant  difference,  the  top  score  (54.01)  for  the 
rhodospin  based  model  was  slightly  higher  than  that  for  the 
/?2AR  based  model.  This  might  come  from  the  weak 
interaction  with  the  EL2  region  in  the  rhodospin  based  model, 
which  was  not  observed  in  the  /32AR  based  one. 

As  in  anibamine,  our  rhodopsin-derived  model  did  not 
indicate  strong  interactions  between  Vicriviroc  and  the  EL2 
of  CCR5,  while  Glu283,  Tyrl08,  and  Ilel98  played  important 
roles  in  the  binding  pocket  (Figure  7G)  in  agreement  with 
recent  site-directed  mutagenesis  results 22f  In  addition, 
Tyr251,  described  as  an  important  residue  in  the  binding 
pocket  of  Vicriviroc  in  the  same  report,  also  had  significant 
interactions  with  the  ligand  in  our  model.  Hydrophobic 
interactions  of  Thr284  with  the  ligand  may  be  significant 
based  on  our  observations  although  this  has  not  been 
examined  experimentally.  In  the  /32AR-derived  CCR5  model 
(Figure  7H),  the  binding  pocket  for  Vicriviroc  shared  some 
features  with  the  rhodopsin-derived  model  with  some  dif¬ 
ferences  arising  from  variations  in  orientation  of  the  amino 
acid  residues  and  the  ligand.  Again,  the  model  indicates 
involvement  of  Tyr251  and  Tyrl08  with  binding,  as  verified 
by  site-directed  mutagenesis  22f  However,  our  observation 
of  Phel09  in  the  binding  pocket  was  contrary  to  the 


mutagenesis  data.  Similarly,  the  roles  of  Gln280,  Trpl90, 
Thr284,  and  Feu255  have  not,  as  of  yet,  been  verified.  While 
the  major  component  in  the  Goldscores  was  the  van  der 
Waals  interactions  between  the  ligand  and  the  receptor,  no 
apparent  difference  was  observed  for  either  model. 

CONCLUSION 

Two  homology  models  of  chemokine  receptor  CCR5  were 
constructed  based  on  the  crystal  structures  of  bovine  rhodop¬ 
sin  and  the  human  /32-adrenergic  receptor,  respectively. 
Anibamine,  the  first  natural  product  CCR5  antagonist,  along 
with  other  three  well-known  CCR5  antagonists  were  docked 
into  these  two  models.  The  binding  pockets  for  these  ligands 
were  characterized  and  compared  with  the  reported  site- 
directed  mutagenesis  studies.  Our  modeling  studies  indicate 
that  the  binding  pockets  for  all  four  CCR5  antagonists  are 
located  in  the  transmembrane  (TM)  part  of  the  protein  near 
Glu283.  This  position  allows  a  putative  ionic  interaction  with 
the  quaternary  nitrogen  that  is  a  common  feature  of  the  CCR5 
antagonists.  While  some  ligands  showed  a  very  strong 
interaction  with  extracellular  loop  2  in  the  rhodopsin 
template-based  model,  others  did  not,  which  may  directly 
influence  their  GOLDScores  (binding  fitness).  In  both  the 
rhodopsin-derived  model  and  the  /32AR-derived  model,  all 
four  ligands  tended  to  adopt  a  similar  binding  locus,  while 
for  each  ligand,  the  binding  modes  and  orientations  were 
somewhat  different  between  these  two  models.  This  is 
certainly  related  to  the  different  3D  conformations  of  the 
TM  domains  and  orientations  of  the  amino  acid  side  chains 
in  the  crystal  structure  template  the  models  were  based  on. 
Most  importantly  for  our  future  work  is  that  anibamine  shares 
many  common  binding  features  with  the  other  known  CCR5 
ligands  in  each  model,  thus  suggesting  that  its  inhibition 
activity  may  stem  from  binding  to  the  CCR5  antagonist 
binding  locus. 

While  from  our  studies,  both  the  crystal  structures  of 
bovine  rhodopsin  and  human  /?2AR  seemed  to  be  reasonable 
templates  to  conduct  homology  modeling  operation  for  CCR5 
receptor,  the  rhodopsin-derived  model  received  more  support 
from  the  reported  site-directed  mutagenesis  studies  for  the 
known  CCR5  antagonists,  while  the  binding  pockets  for  those 
ligands  identified  from  the  /32AR-derived  model  need  to  be 
further  verified  by  more  extensive  experimental  studies. 

Anibamine  is  an  interesting  lead  compound  for  anti-HIV 
drug  discovery.  It  is  the  first  and  only  natural  product  CCR5 
antagonist  identified  such  far,  and  because  of  its  unique 
skeleton,  the  recently  reported  flexible  total  synthesis  of 
anibamine  opens  up  a  number  of  research  directions.  In 
particular,  the  pursuing  of  the  next  generation  molecular 
design  based  on  the  unique  binding  features  of  anibamine 
may  yield  a  clinically  relevant  CCR5  antagonist.  The  binding 
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mode  study  of  anibamine  in  CCR5  homology  models 
reported  herein  will  be  beneficial  to  the  future  study  in  novel 
drug  design  and  development. 

COMPUTATIONAL  METHODS 

The  sequence  of  human  chemokine  receptor  CCR5  was 
obtained  from  http://www.expasy.org  (Swiss-Prot  and  TrEM- 
BL),  and  the  crystal  structures  of  bovine  rhodopsin  and 
human  /32AR  were  downloaded  from  http://www.rcsb.org. 
All  computations  were  performed  on  SGI  Octane2  worksta¬ 
tions  (Silicon  Graphics,  Inc.,  Mountain  View,  CA). 

Homology  modeling  was  carried  out  using  the  Homology 
and  Discover  modules  of  the  Insightll  package  (Accelrys, 
San  Diego,  CA).  The  sequence  alignment  was  conducted 
automatically  by  applying  different  scoring  matrices,  gap 
penalty,  and  gap  length  penalty  parameters  in  order  to 
achieve  the  highest  sequence  homology.  Then  the  coordinates 
for  the  amino  acid  residues  in  the  structurally  conserved 
regions  of  CCR5  were  assigned  following  the  corresponding 
residues  in  the  CCR5  receptor.  The  coordinates  for  the  amino 
acid  residues  of  the  CCR5  receptor  in  the  structurally 
nonconserved  regions  were  obtained  by  searching  the  protein 
database  with  the  InsightII/Homology/Loop_Search  module. 
After  the  completion  of  the  amino  acid  residue  coordinates 
assignment,  the  conformation  of  the  whole  protein  was 
minimized  with  the  backbone  constrained.  A  brief  minimiza¬ 
tion  (2,000  iterations)  and  dynamics  simulation  (10,000  step, 
1  fs  each  step)  was  conducted  to  relax  major  steric  clashes 
introduced  during  the  coordinate  assignment.  Further  vigor¬ 
ous  simulation  (minimization  10,000  iterations  and  dynamics 
100,000  steps,  100  ps  total)  was  conducted  to  produce  an 
energetically  reasonable  conformation.  The  CFF91  force  held 
was  applied,  while  the  minimization  was  terminated  until 
the  maximum  derivative  is  less  than  0.001  kcal/mol. 

All  the  small  (ligand)  molecules  were  built  using  Insightll/ 
Builder  Module.  Minimizations  with  steepest  descent  fol¬ 
lowed  by  conjugate  gradient  were  performed  to  generate  the 
lowest  energy  conformation  for  each  ligand  studied  (CFF91 
force  held  and  default  termination  values  were  adopted). 
Then  a  molecular  dynamics  simulation  was  performed  (an 
equilibration  phase  of  1000  fs  at  300  K,  followed  by  a 
collection  phase  of  5000  fs  at  the  same  temperature)  to 
further  study  the  small  molecule  conformation.  We  used  the 
lowest  energy  conformation  of  the  small  molecules  in  water 
extracted  from  the  5  ps  molecular  dynamics  calculations,  as 
the  initial  conhguration  for  docking  into  the  proposed  binding 
site  of  the  CCR5  receptor.  The  ligand  was  modeled  in  its 
nitrogen-protonated  form  when  applicable. 

The  interactive  docking  study  was  conducted  in  the 
Discover  module  of  Insightll.  The  ligands  were  docked  in 
the  upper  portion  of  transmembrane  region  of  the  receptor, 
following  orientations  reported  in  the  literature.22  Each 
ligand— receptor  complex  was  minimized  in  the  gas  phase 
hrst  with  the  backbone  of  the  receptor  fixed,  but  all  the  side 
chain  atoms  were  left  unconstrained.  The  optimized  confor¬ 
mation  was  then  used  as  the  initial  conhguration  for 
molecular  dynamics  simulations.  For  each  complex,  a  short¬ 
term  steepest  descent  energy  minimization  (500  iterations) 
and  a  dynamics  simulation  (1000  step,  1  fs  each  step)  was 
conducted  to  validate  the  initial  docking  conformation.  This 
was  followed  by  another  simulation  minimization  (5000 


iterations),  and  100,000  steps  of  dynamics  were  conducted 
with  an  initial  2000  steps  equilibration.  The  total  simulation 
time  was  102  ps.  Again  the  CFF91  force  held  and  default 
termination  value  (0.01  kcal/mol)  was  adopted. 

Automated  docking  was  conducted  with  GOLD.55  Default 
parameters  were  adopted  except  that  the  docking  locus  was 
dehned  as  10  A  around  the  OE1  of  Glu283  on  helix  7  of 
CCR5.  Based  on  the  htness  scores  and  binding  orientation 
of  the  ligands,  the  best  solution  of  GOLD  binding  (among 
top  ten  solutions  that  were  recorded)  was  selected  and  merged 
into  the  structure  of  CCR5  receptor.  The  combined  structures 
of  the  ligand  and  CCR5  receptor  were  minimized  together 
in  order  to  allowing  the  receptor  to  adopt  the  presence  of 
ligand.  After  that,  the  ligand  was  removed  from  the  complex, 
and  the  resulting  structure  was  evaluated  by  PROCHECK 
again.  The  docking  results  were  summarized  in  Table  3. 
Finally,  all  the  models  were  optimized  by  molecular  dynam¬ 
ics  simulation,  using  SYBYL7.3  with  default  parameters 
except  for  the  following:  length  -  100,000  fs,  snapshot  every 
25  fs;  time  step  -  1  fs,  dielectric  constant  4.0. 

The  PDB  files  of  the  homology  models  are  available  from 
the  corresponding  author  upon  request. 
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